Novel Phenomena in Charged Bose Liquid 
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We investigate charged Bose liquid immersed in uniform background charge at zero temperature. 
Novel phenomena, such as oscillatory shielding of external localized electric charge, rotons and charge 
density waves (charge stripes in two dimensions), occur in any dimensions. Oscillatory shielding is 
caused by mixing between scalar boson exchange and Coulomb interactions, which mediate opposite 
forces. On the other hand, rotons and charge density waves are due to attractive local self-interaction 
of bosons. Rotons can be regarded as a finite size charge density wave packet without any back flow. 
We also comment on charge stripes observed recently in cuprates and nickelates. 
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Recently there have been several experimental obser- 
vations of charge and spin stripes in high temperature 
superconductor materials Jl|-§| . There have been several 
theoretical discussions about the origin and implications 
of these stripes [Q. Experiments show on-set of charge 
order before spin order, showing that charge stripe is pri- 
mary and spin stripe is secondary. The superconducting 
transition temperature is either zero in compounds with 
a static wave of charge density, or significantly reduced 
when dynamical correlations are observed. It would be 
interesting to find a simple model which produces charge 
stripes and dynamical correlations. 

In this paper we investigate charged Bose liquid im- 
mersed in uniform external electric charge at zero tem- 
perature. We describe charged Bose liquid by a complex 
scalar field with nonlinear local self-coupling and gauge 
coupling with photons. By analyzing linear fluctuations 
around homogeneous ground state, we show that novel 
collective modes, like oscillatory shielding of external lo- 
calized charges, rotons and charge density waves, occur 
in any dimensions. We identify physical reason behind 
these phenomena. Our model can be regarded as a sim- 
ple toy model which provides such rich phenomena, and 
may also shed some light on current experiments. 

Since unique is the rest frame of uniform external 
charge, there is no Lorentz or Galilean symmetry in our 
theory. There have been many works on this system with 
attention on vortex dynamics K|- Vortex dynamics is 
particularly simplified in the self-dual model 0, where 
magnetic flux vortices are shown not to carry any intrin- 
sic angular momentum. 

The presence of uniform background charge leads to 
nontrivial consequences. The average charge density of 
bosons is determined by the condition of electric charge 
neutrality rather than by the chemical and potential en- 
ergy. For this reason, the (short-range) self-interaction 
of bosons does not have to be repulsive for the system 
to be stable. If this interaction is weakly attractive, 
low-energy excitations from a uniform ground state are 
fluctuations of charge density with a characteristic wave- 
length. These excitations can be called rotons in analogy 
with those in superfluid helium. As the strength of at- 
traction increases, the frequency of this mode vanishes 
and the uniform ground state becomes unstable, form- 



ing charge density waves. In superfluid helium nonlocal 
attractive self-interaction of helium atom plays a crucial 
role in the existence of rotons. (For a review of cur- 
rent understanding of rotons in superfluid helium, sec 
Rcf. @.) 

A closely related question is the screening of a probe 
charge. It is normally expected that a localized external 
charge is fully shielded by charged bosons; the amount of 
total charge inside a sphere with the external charge in its 
center falls off exponentially with the radius. However, in 
the case we study, the shielding can be oscillatory. This 
originates from mixing of two opposite forces mediated 
by scalar boson exchange and Coulomb interactions. 

We find a link between oscillatory shielding and the 
existence of a roton dip in the excitation spectrum. As 
stressed by Feynman || , the roton minimum is associated 
with a short-range crystal-like structure in liquid helium. 
In that case, the density-density correlation (p(v)p(O)) 
is a non-monotonic function of the distance, an effect 
similar to the oscillatory screening of charge. 

Within the scope of our model, a roton wave packet 
can be visualized as a localized charge density wave. As 
discussed above, it originates from the attractive self- 
interaction of bosons. While the wave packet itself is 
at rest (zero group velocity), charge density waves are 
moving at the roton phase velocity. Curiously, there is 
no net flow of charge associated with this motion and 
therefore it is unnecessary to attach a vortex ring to this 
wave packet in order to produce a back flow. The idea 
of back flow goes back to Feynman and Cohen || , who 
tried to modify Feynman's original ansatz to enforce the 
conservation of current. In our theory the charge current 
is exactly conserved from the start and so there is no 
need for back flow. At least our result is consistent with 
the current view of rotons in superfluid helium . 

We describe charged Bose liquid by a bosonic scalar 
field yb with a local self-interaction and the electromag- 
netic coupling. This liquid is immersed in a uniform back- 
ground charge p e > 0. The Lagrangian for charged Bose 
liquid is 

C = -^-F u „F^ + ii/,*(do-iqA )il> (1) 
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where h = c = 1. The Gauss law constraint is 

d t F 0l + 4Tcq(p - Pe ) = 0. (3) 

As the charged background is static, homogeneous and 
uniform, there are conserved energy and linear and angu- 
lar momenta. The Noether angular momentum density 
is not gauge invariant and needs a modification, which 
plays a crucial role in vortex dynamics |5) . 

The homogeneous configuration satisfying the field 
equations is given as 

(A M ) = and (iP) = Vple- Mt , (4) 

where a = dU/dp\ Pe . Total charge density is equal to 
zero and the system is charge neutral. We approach this 
problem by the semiclassical approximation around this 
uniform background. To find out elementary excitations 
on this homogeneous configuration, we study linear fluc- 
tuation equations in the unitary gauge Arg ip = 0, and so 
(Aq) = a/q. In the Fourier mode SipjSA^ ~ e -*wt+i P l z^ 
the linearized field equations become 

2 

-(|- + 2Ap e )(5V + qy/pl&Aa = 0, 
Zm 

p 2 5A + ujp^Ai + %-Kq s /pl5ij) = 0, 

(u 2 - p 2 - m 2 A )5A t + p^SAj + up^Ao = 0, (5) 

where A = d 2 U/dp 2 \ Ps and m 2 A — Anp e q 2 /in. The pa- 
rameter A is the strength of the local self-coupling for 
linear fluctuations. 

From these equations, we get the dispersion relation 
for the transverse modes, 

(uj 2 -p 2 -m 2 A )e ljkP iA k =0. (6) 

This Anderson-Higgs mechanism explains the Meissner 
effect with the London penetration length Xl = l/ m ^- 
The dispersion relation for the mixed mode of Sip, 5Aq 
and p l 8Ai is 

For the time-independent fluctuations, we get the static 
relation 

p 2 = -2m\p e ± 2m^(Xp e ) 2 - m\. (8) 

This describes how the system responds to any time- 
independent external source. 

Depending on the range where the parameters of the 
theory lies, we will see that several collective modes are 
allowed by Eqs. (^) and (^), which are summerized as 
follows: 

Xp e > ttia'-i exponential shielding, (9) 
\Xp e \ < ttia; oscillatory shielding, (10) 
> Xp e > —rnA] rotons, (11) 
Xp e < —iriA', charge density waves, (12) 
Xp e « —rriA] Wigner crystal. (13) 



When the electromagnetic coupling constant e 2 van- 
ishes, this theory has a smooth limit. It describes a neu- 
tral Bose liquid with local self-interaction. We know well 
that with repulsive self-interaction A > 0, the homoge- 
neous configuration is stable with positive pressure. El- 
ementary excitations in that case are massless phonons 
with sound speed v 2 — Xp e /m. If we have introduced 
a proper nonlocal self-interaction as in superfluid he- 
lium, the roton mode would have appeared. There is 
no massive mode but there exists a natural length scale, 
i.e., the coherence length £ = l/y/4mXp e . With attrac- 
tive self-interaction A < 0, the configuration is unsta- 
ble. Bosons aggregate and form droplets, or Q-balls, 
with mean charge density p Cl at which U(p) takes the 
minimum value. 

Let us first study how the system responds to the static 
localized probe charge. Far away from this source, the 
field configuration deviates little bit from the uniform 
configuration (0). Since the charge couples to A$, Eq. (||) 
would govern the behavior of such a deviation. There are 
three interesting ranges of the parameters from Eq. (|§). 
When the coupling is strongly repulsive so that the pa- 
rameters lies in the range (||) , the solution of Eq. (||) is 
pure imaginary. Thus the local source is shielded expo- 
nentially. However, the solution of Eq. (||) becomes com- 
plex when the parameter lies in the range (|l0|) . Thus the 
fields would fall off exponentially but the sign of the field 
would oscillate, like the amplitude of underdamped har- 
monic oscillators. We will call this as oscillatory shield- 
ing. 

To find the physical origin of such oscillatory screen- 
ing, let us recall the old lore. For two like sources in 
large separation, the virtual scalar boson exchange leads 
to the attractive Yukawa force, but the virtual photon 
exchange leads to the repulsive Coulomb force. For the 
opposite charges, the scalar force is repulsive and the 
Coulomb force is attractive. This can be seen directly by 
obtaining the interacting potentials from the correspond- 
ing Lagrangians, 

C 1 = -^(dJ) 2 -^m 2 J 2 +p(x)f, (14) 
C 2 = ^(d l A Q ) 2 + ^m 2 A 2 +p(x)A , (15) 

where p(x) — qi5 3 (x — xi) + q2$(x — xi). Here we have 
neglected the time dependence and scaled the field vari- 
ables. The sign of the field propagator is also associated 
with the sign of the exchange force. Thus the A^ ex- 
change would lead to the similar force as the scalar boson 
exchange. We already know this is true as two parallel 
currents are attractive and two opposite currents are re- 
pulsive. (The magnetic flux of vortices is not a direct 
source of Ai unlike the current. The force between fluxes 
by the Ai exchange is more like the Coulomb case.) When 
there is a mixing between Aq and / by a Lagrangian 
£ 3 = p 2 fA , the corresponding static relation of the to- 
tal Lagrangian C\ + £2 + £3 becomes 

(p 2 + m 2 )(j5 2 + m 2 )+/i 4 = 0. (16) 
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When 2p? > \m\ — m 2 ,] so that the mixing is strong 
enough, the solution p 2 changes from negative real num- 
bers to complex number. This is how oscillatory shielding 
appears. 

More interestingly, when the parameters lie in the 
range (p^), the solution of Eq. (JsJ> becomes real num- 
bers. The fluctuations would not die out. Especially near 
Xp e « -mj, the charge density would be oscillating with 
a characteristic wave length 

2ir , 

ICDW = /f) (I'J 

This describes charge density waves. (In two dimensions, 
it would be charge stripes.) In our case m 2 and m\ of 
Eq. ( pif ) can be negative and so the charge density wave 
mode is possible. 

Let us now consider the time-dependent fluctuations 
described by Eq. (fjj). The slope of u)(p) at p = becomes 
negative when the self-interaction is attractive so that the 
parameters lie in the range (pd[). It has the minimum at 

p, = ^2m(-Xp e ), (18) 

with the value 

= jm\ - X 2 p 2 e . (19) 

As long as Xp e > — tua, the mode at the minimum has 
the positive energy and is called rotons. 

For a roton mode, we have to choose a particular direc- 
tion, p*. We can form the roton wave packet by super- 
posing the plane waves with an envelop function T . Since 
the dispersion relation distinguishes the parallel and or- 
thogonal momenta with respect to p*, T can depend on 
parallel and orthogonal momenta independently. Espe- 
cially if T is peaked at p„, the group velocity would be 
zero. 

Now one may wonder what happens to charge density 
and current in the roton wave packet, especially whether 
there exists any back flow. We first note that our theory 
works in any spatial dimension. Thus we first consider 
the system in one dimensional space. When the wave 
packet size is large, the major contribution to charge den- 
sity and current comes from the p* plane wave. Inside the 
wave packet, both positive and negative charge regions 
of wave length A* = 27r/p* moves forward with phase ve- 
locity v* = uj*/p*. Thus, there is no net charge transfer 
from one end to another end of the wave packet, and no 
need for back flow of charge. This is consistent with the 
fact that the current conservation holds for linear fluctu- 
ations. The charge density wave inside the roton wave 
packet is exactly like sound wave in fluid, which does not 
transport any mass density and needs no back flow. 

Coming back to the two or three dimensional prob- 
lem, we note that the envelop function can be chosen 
to be Gaussian along the parallel and orthogonal direc- 
tions, independently. Then the above one-dimensional 
argument shows that there is no back flow even in these 



cases. If there is any back flow, it may manifest itself as 
vortex ring, or smoke ring around the roton wave packet, 
which is clearly nonperturbative and contradicts the per- 
turbative nature of our fluctuations. Smoke ring is also 
impossible in one-dimensional space. 

If the interaction between bosons become more attrac- 
tive so that Eq. (|l2| ) holds, the roton mode becomes un- 
stable and charge density waves appear. The homoge- 
neous field configuration would be no longer the ground 
state configuration. Since the charge density wave mode 
is characterized by a single length scale Icdw, the ground 
state is made of a static plane wave of the given wave 
length, pointing a particular direction, say z. This low- 
ers the gradient energy. Thus the translation symmetry 
is broken along the z axis but is invariant along the x-y 
plane. The rotational symmetry is also partially broken 
to the z directional rotation. There would be a Gold- 
stone boson mode which modulate the position of the 
crest of charge density waves. The amplitude of charge 
density waves would be fixed by the nonlinear interac- 
tion. The exact semiclassical description of the ground 
state needs the evaluation of the full nonlinear field equa- 
tions, which will tell us whether the background charge 
is exposed at the bottom of the wave. (This may tell us 
whether this phase is superconducting or metallic along 
the z-direction). 

When the attractive interaction becomes still stronger, 
the unstable modes can have a range of momentum. We 
have to treat the full nonlinear field equation to find the 
exact ground state field configuration. However, one can 
guess the ground state configuration in this case by start- 
ing from the other end where the electromagnetic cou- 
pling is negligible. If we turn off the electromagnetic cou- 
pling, charged bosons will aggregate to many droplets, 
or Q-balls, and then they will grow. When we turn on 
the electromagnetic coupling, they cannot grow forever. 
Rather Q-balls of suitable size would arrange themselves 
into a crystal structure to lower the Coulomb energy. 
This configuration is analogous to the Wigner crystal of 
electrons in a uniform positive charged background. In 
this case the translation and rotational symmetries are 
broken to a discrete crystal group. 

The recently discovered charge stripes and dynamical 
correlations can be regarded as charge density waves and 
rotons in our context. While their dynamical origin may 
be far complicated, charged Bose liquid considered here 
may capture some physics (aside from spin and lattice 
effects). As the dopping increases in these materials, the 
carrier charge density changes. There are several possi- 
bilities of behavior. Bosons can be repulsive or attractive 
in small density, though they are expected to be repulsive 
in high density. 

At very low density itia ~ *J~pl 3* | Ap e j and therefore 
the charge shielding is oscillatory with the wave vector 
p rs (i ± l)y/mmA. Experimentally it may be difficult 
to tell this type of shielding from purely exponential be- 
cause the imaginary part of p equals the real part. If A 
is positive and approximately constant, the exponential 
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screening appears in high density. 

Second possibility is that A is negative in small den- 
sity and so the roton mode presents in low density. As 
the density p e increases, the roton dip becomes more 
prominent. Being lowest energy excitations of the sys- 
tem, rotons deplete the condensate density at tempera- 
tures equal to or above their energy, like dynamical cor- 
relations in real material. In higher density the repulsive 
interaction turns on and the roton mode disappears. 

If the attractive interaction remains to high enough 
densities, the roton mode becomes unstable at some crit- 
ical density and the system may not be superfluid any 
longer. At a higher density, the reduced attraction can 
no longer maintain a static charge density wave and su- 
perfluidity comes back. At the moment of instability the 
wavelength of charge density waves is determined by the 
boson mass m and the photon mass ttia (the London 
penetration depth). If we take m equal to two electron 
masses and the London length Al = 2000A, a typical 
experimental value, the resulting CDW wavelength (|lT 
is of order 
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8.5k. 



(20) 



This is of the same order as the experimental value ~ 
15 A, which is encouraging. (Note that this analysis is 
bases on linear fluctuations and so the Meissner effect 
is still true. But this may not be true once the exact 
semiclassical ground state is found.) 

It is surprising that our simple system shows many fas- 
cinating phenomena, as oscillatory shielding, rotons and 
charge density waves. Our work suggests that these phe- 
nomena are quite generic. We now have tools to explore 
quantum mechanics of these collective modes. We can, 
for example, calculate spin of rotons in our model. It 
would be also interesting to analyze the system by us- 
ing the method of many body physics. It is known that 
when A = the homogeneous configuration changes to 
the Wigner crystal in high enough density due to the 
Coulomb repulsion |uj, which we do not see in our ap- 
proach. 

Even though we have not considered here, there exist 
also magnetic flux vortices in the theory. The naive di- 
vision between type I and type II should be redrawn by 
studying the long-distance interaction between magnetic 
flux vortices in the exponential shielding region. 

It would be most interesting to find out whether our 
consideration has any direct bearing on real superconduc- 
tors, whose physics was investigated in context of charged 
Bose liquid jll] . Recently there has been some proposal 
of more realistic field theoretic models of the BCS super- 
conductors at zero temperature It would be inter- 
esting to find out whether our consideration applies to 
that case. There are also numerous works relating super- 
conductivity and charge density wave in one-dimensional 
systems We hope that our work can provide some 
useful insights here as well. 
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